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Abstract — In this paper, convolutional network coding 
is formulated by means of matrix power series representa- 
tion of the local encoding kernel (LEK) matrices and global 
encoding kernel (GEK) matrices to establish its theoretical 
fundamentals for practical implementations. From the en- 
coding perspective, the GEKs of a convolutional network 
code (CNC) are shown to be uniquely determined by its 
LEK matrix K{z) it Kq, the constant coefficient matrix 
of K{z), is nilpotent. This will simplify the CNC design 
because a nilpotent Kq suffices to guarantee a unique set 
of GEKs. Besides, the relation between coding topology 
and K{z) is also discussed. From the decoding perspec- 
tive, the main theme is to justify that the first L + 1 terms 
of the GEK matrix F{z) at a sink r suffice to check whether 
the code is decodable at ? with delay L and to start decod- 
ing if so. The concomitant decoding scheme avoids dealing 
with F{z), which may contain infinite terms, as a whole 
and hence reduces the complexity of decodability check. It 
potentially makes CNCs applicable to wireless networks. 

I. Introduction 

Network coding was formally introduced by fl\ . Later, lin- 
ear network coding was proved to be able to achieve the opti- 
mal data transmission rate in an acyclic multicast network |2|, 
and an algebraic approach to linear network coding was pre- 
sented in |[3|. Since then, a rich literature on Unear network 
coding has emerged, and a wide variety of applications have 
been developed. 

Over a network with cycles, the propagation and encod- 
ing of sequential data symbols naturally convolve together; 
the propagation delay thus becomes an essential factor in net- 
work coding. To cope with cyclic data transmission, convo- 
lutional network coding was introduced in |T,3''6"'T0l . It is a 
form of linear network coding which deals with the pipeline 
of messages as a whole rather than individually. Under the 
assumption of unit-delay edge transmission and local encod- 
ing kernels chosen from a finite field, an algebraic framework 
for a CNC was formulated in This framework is general- 
ized in |6| for those CNCs in which along every cycle there 
is at least one delay, and the LEKs are chosen among rational 
power series over the symbol field. The ring of rational power 
series was justified therein to be the proper algebraic structure 
for a CNC. More recently, the framework was mathematically 
extended in |jTT| for such CNCs in which LEKs determine a 
unique set of GEKs. In each framework, the main theorem 
guarantees the existence of an optimal CNC under the respec- 



tive assumptions such that the GEK matrix at every receiver 
has full rank. In order to efficiently construct such an optimal 
CNC over a cyclic network, a polynomial-time algorithm was 
proposed in |9|. On the other hand, a method was introduced 



in 1 1 1 1 to adapt any acyclic algorithm for the construction of 
an optimal CNC on a cyclic network. 

Previous studies on convolutional network coding examines 
every data unit of time-multiplexed data symbols as a whole, 
and manipulates the LEK matrix over an algebraic structure 
of data units. In this manner, a CNC is nothing but a linear 
network code over this algebraic structure. Thus, the classical 
field-based algebraic framework in f3^ for acyclic networks 
could be applied for a CNC with some additional assumptions 
and modifications. On the other hand, this facility obscures 
the implementation aspects of CNCs. 

On the encoding side, a CNC is generally deployed in every 
encoding node by its LEKs. For acyclic networks, GEKs can 
be uniquely determined by LEKs f5\, and the practical fea- 
sibility is also assured. However, over a cyclic network, the 
GEKs may not be uniquely deduced from the LEKs (see ex- 
amples 3.2 and 3.3 in |5|). Sometimes, the code may not be 
feasible even if the GEKs can be uniquely deduced because 
the deduced GEKs do not satisfy the rational form. In this 
paper, we firstly start the study of a practically feasible CNC 
from a new approach, that is, the power series representation 
of a matrix. 

We adopt the ring F[(z)] of rational power series over the 
symbol field F as the ensemble of data units for a CNC, where 
z is the unit time delay. Just as every rational power series over 
F can be written as J2t>o ^tz*, kt G F, we shall represent a 
matrix K{z) over F[(z)] as a rational power series J2t>o ^tz*, 
where Kt is the matrix over F. For example. 
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With such a novel representation, we can characterize the 
field-based conditions for a CNC 

• to determine a unique set of GEKs, which is a prerequi- 
site for data propagation; and 

• to be practically feasible. 

The related work on the encoding side of CNCs is summa- 
rized in Fig{T| and will be explored in Section III, after some 



related fundamentals on CNCs are reviewed in Sec. II. In de- 
tail, we first consider the implementation of a CNC and find 
that a CNC is practically feasible if and only if the operations 
at intermediate nodes are finite and causal, which is also re- 
ferred in |6 |. In order to characterize the property of finite and 
causal operations, we give the definition of expandability of a 
matrix. Further more, logical inconsistency exists if there is no 
partial encoding order along the cycle. So we define an Encod- 
ing Topology (ET) w.rt. Kq to illustrate the encoding order. 
As a result, it is shown that the GEKs can be uniquely formu- 
lated from LEKs when the ET w.r.t. Kq is acyclic, in which 
case Kq is nilpotent and results in expandable GEKs. This im- 
plies whether a CNC is practically feasible is only determined 
by Kq, which simplifies the CNC encoding design. Finally, 
we present some equivalent conditions in the same section. 
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sion of data symbols. Then we provide several sufficient and 
necessary conditions for the code to be decodable at this re- 
ceiver with delay L, which only involve coefficient matrices 
Fo , Fi , • • • , over F. These conditions are field-based and 
hence greatly reduce the computational complexity in code de- 
sign. The relations among different conditions for decodabil- 
ity of a CNC are summarized in Figj2j 



A prdtitally feasible CNC with GEK 
nfiatrix Fl{z) for sini< r is decodable 
at r with delay L. 
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There is a matrix Diz) over F[(z)] 
such that F(z)D{z) = zH 
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mnkiFi) — rank(Fi-i} = j: 



Figure 2: the relation among decoding conditions of a CNC. 
Here Xq{z) is the sequence of source symbol vectors gener- 
ated from time unit to i and j/q (z) the sequence of received 
symbol vectors at receiver r over the same period. Based on 
matrix power series representation, F^{z) is the first L + 1 
terms of F{z) and is defined by (mi. 



Figure 1: the organization and relations among the conditions 
to uniquely determine GEKs based on the LEK matrix K{z) = 

On the decoding side, in previous work, a CNC is decod- 
able at a receiver if and only if the GEK matrix over F[(z)] at 
this receiver has full rank. Moreover, in order to decode with 
delay L, the current method (See |5| for example) is to com- 
pute a matrix D{z) over F[(z)] such that F{z)D{z) — z^I, 
where F{z) is the GEK matrix at this receiver and D{z) is 
called the decoding matrix. However, due to the existence of 
cycles in the network, GEKs may involve infinite terms. In the 
decentralized deployment of a CNC, it is infeasible for a re- 
ceiver to judge whether the code is decodable after collecting 
all information on its GEK matrix. Even if the GEK matrix 
is known by the receiver as in the centralized deployment, the 
computational complexity to check full rank of F{z) will be 
very high. 

In section IV, based on matrix power series representation 
of the GEK matrix F{z) for a receiver and motivated by p3] , 
we give another definition of decodability at a receiver with 
delay L, which captures the feature of sequential transmis- 



II. Definitions and related works 
A Notation 

A communication network is modeled as a finite directed 
graph Q = (V, £) with possible cycles. A directed edge rep- 
resents a noiseless communication channel transmitting a data 
symbol per unit time. Assume there is a unique source node, 
denoted by s, in the network. The symbol alphabet is a finite 
field denoted by F. In each unit time t, the source s generates 
a message, which consists of a fixed number lj of symbols and 
is presented by an w-dim row vector xt over F. For every node 
V, denote the set of incoming channels by In{v), and the set of 
its outgoing channels by Out{v). For technical convenience, 
we assume that In{s) consists of lu imaginary channels. An 
ordered pair (d, e) of channels is called an adjacent pair when 
there exists a node v with d E Iniv) and e e Out{v). 

Via a CNC, a stream of messages propagates from s through 
a time-invariant convolutional encoder at every node, and it is 
represented by an oj-dim row vector of x{z) — X]t>o ^tz*", 
where Xt G F" and z is the time variable. Every entry in x{z) 
belongs to the principal ideal domain (PID) F[[z]] of power 
series. Such power series in F[[z]], which can be written in 



the form of p{z)/{\ + zq{z)), where p{z) and q{z) are poly- 
nomials, are called rational power series. Denote the PID of 
all rational power series by F[(z)]. A square matrix A{z) over 
F[(2;)] is invertible iff there exists a matrix C{z) over F[(2:)] 
such that A{z)C{z) = /. The determinant of A{z) can be 
indicated by A&i(A{z)) = oq + aiz + 022;^ + • • • G F[(z)]. 
A sufficient condition for A{z) to be invertible over F[(z)] 
is the nonzero value of oq. Denote the adjoint matrix of 
A{z) hy A* {z). Thsn, A{z)A*{z) = det(A(z))/. If oq is 
nonzero, we have A{z)^^ = A* (z) / det{A{z)) and it can 
be expressed as a positive power series of z. Otherwise, it 
is not clear for us whether A{z) is invertible. Assume that 
det{A{z)) = z\at + at+iz + •••), and at 7^ 0, < > 0, the 
invertibility of A(z) depends on whether the entries of A*{z) 
have the common factor z* or not. 

We next define the matrix power series representation of a 
matrix over F[(z)]. 

Definition 1. For an m x n matrix A(z) over F[(z)], its ma- 
trix power series representation is X^t^o ^t-^*' where At is an 
m X n matrix over F in which the {i,jY^ entry is equal to the 
coefficient of term z* in the {i,jY^ entry in A{z). Two matrix 
power series are equal if and only if the coefficients of z* are 
equal for any t. 

Definition 2. Let A{z) be a square matrix over F[(z)], and 
B{z) a function of A{z) in the form ofY^'^^ CtA*{z), where 
Ct,t > 0, are square matrices over F. We say that B{z) = 
X]t>o ^tz* is expandable if for all t, the matrix coefficient 
Bt can be written as the sum of finite terms in the form of 

Definition [2] is motivated by the characterization of prac- 
tical transmission operations. That is, all the operations at 
each intermediate node and sink are finite and causal for each 
time slot t. As will be justified in Theoremj2] if B{z) = 
X]t>o^('^)*' ^^^^ B{z) is expandable if and only if Aq is 
nilpotent. 

B Fundamentals on CNC 

In this section, we review the local and global encoding kernel 
descriptions of a CNC. 

Definition 3. (Local description) An uj-dim F-CNC on a com- 
munication network with possible cycles consists of an element 
kd,e{z) € IP'il^)]. called the local encoding kernel (LEK), for 
every adjacent pair {d, e). 

Definition 4. (Global description) A set of global encoding 
kernels (GEKs) for an uj-dim F-CNC with LEKs kd^e{z) is 
an assignment of an uj-dimensional column vector fe{z) over 
V[{z)]for every channel e such that: 

1) fe{z) = J2dein{v) kdAz).fd{z) when e E Out{v). 

2) The vectors fe{z), e G In(s), form the natural basis of 
the free module F[(z)]". 



Over an acyclic network, it is equivalent to define a CNC on 
GEKs and LEKs. However, it is not the case in a cyclic net- 
work |5|. Sometimes LEKs determine multiple sets of GEKs 
whereas it is also possible for different LEKs to yield a same 
set of GEKs (See Fig. 5 and Fig. 7 in Sec. Ill for example). 
Normality of a CNC was introduced in 1 1 1 1 . 

Definition 5. A CNC is said to be normal iff LEKs determine 
a unique set of GEKs. 

The following lemma justifies normality of a CNC as a pre- 
requisite for data transmission via the code. 

Lemma 1. For a normal CNC, the symbol ye,t trans- 
mitted over each channel e G Out{v) at time t is 

J2deln{v) (j^T^a kd.e,TVd,t—r^- 

Proof: Let source s generate a message x{z) — ^ xtz*, 

t>o 

which is an w-dimensional row vector over F[(z)]. Through 
a CNC, each channel e carries the power series x{z) ■ fe{z) of 
data symbols. That is, 

ye{z) = x{z) ■ fe{z) 

= x{z) ■ Y kd,e{z)fd{z) 
deln(r) 

kd,e{z){x{z) ■ fd{z)) 

deln(r) 

= Y kdAz)yd{z) 

d&In(r) 

In matrix power series representation, T/e(z) = X]t>o Ve,tz^ 
and kd^z) = Et>o kd,e,tz*. Then 

0<T<t 

and 

ye,t = ^ ^ kd,e,Tye,t-T | • 
deln{r) \t=0 J 

□ 

Let n be the number of channels in the network. For a nor- 
mal CNC, denote by K{z) the nxn matrix [fcd,e(-z)]d,eG£> ™d 
F{z) the u! X n matrix [fe{z)]e££- Let Hs represent the uj x n 
matrix [I^j 0]. Then we have the classical equation from |3 1 

F{z) = H, + Fiz) ■ K{z) 

and we have 

F(z){I,,~K{z))^Hs (1) 

According to ([T]i, if /„ — K{z) is invertible, then the code 
is normal. A sufficient condition for the invertibility of /„ — 



K{z) is that / + K{z) + K'^{z) -\ is expandable, i.e., 

(/„-if(z))(^i?,z*) = {I-Kiz)){f2K\z)) 

t t^O 

oo oo 

= {Y,K\z))-J2k\z) 

oo oo 

= I + J2k\z)-J2k\z) 

t=l t=l 

= / 

where J^t ^t-^^* = St^o K'^(z). That is, /„ — K{z) has an 
inverse B{z) = J^t ^tz*- So we have the following proposi- 
tion. 

Proposition 1. /„ — K{z) is invertible if I + K{z) + (z) + 
■ ■ ■ is expandable. 

However, /„ — K{z) is still possible to be invertible when 
/ + K{z) + K'^{z) + • • • is not expandable. An example is 
shown in Fig. [9]in Section III. 

C Practical feasibility of a CNC and Encoding 
topology 

Normality of a CNC is not a sufficient condition for practi- 
cal implementation of the code because we cannot wait for the 
whole GEKs to start decoding at the receiver, especially in the 
case of infinite GEKs. In ]6p I), the causality of a CNC is fur- 
ther defined to guarantee physical implementation. However, 
in decentralized practical application, receivers can only get 
the GEKs term by term. According to this term by term trans- 
formation feature, below we formally justify that causality is a 
necessary and sufficient condition for a code to be practically 
feasible in terms of encoding topology (ET). 

Definition 6. A normal ¥-CNC is said to be practically fea- 
sible iff each node v can calculate the symbols to be trans- 
mitted over its outgoing channels at time slot t by its previ- 
ously received symbols and LEKs kd,e.T, where d € In{v), 
e G Out[v), and t < t. 

Definition |6] requires that every operation at intermediate 
nodes is causal and finite. Similar to linear network codes, we 
would like to find a partial order to define the causal encod- 
ing operation. For every CNC on a network, we define ET 
with respect to the LEK matrix K{z) to be a directed graph 
in which £ is the node set and there is a channel from d to e 
iff kd.e(z) 7^ and take kd.e{z) as the coding weight over the 
link. 

Example 1. Assume that we multicast two messages from 
source S to both the nodes X and Y with given LEKs depicted 
in Fig. |4] Easily, we have the ET w.r.t. the LEK matrix K{z) 
shown in Fig. |3]and the ET w.rt. in Fig. |5] 

To perform a CNC on a network we need an encoding order 
-<i on the node set for all time slots. We call a set of encoding 
orders time invariant if the encoding orders at all time slots 




Figure 3: A convolutional network code with Kq nilpotent 




Figure 4; The ET w.r.t. K{z) in Fig. |3]is acyclic. 



are the same. A CNC is practically feasible if it is practically 
feasible w.rt. a set of encoding orders. ET implies the coding 
relations among edges, and we can designate a partial order 
when the ET is acyclic. Further, we find an acyclic ET w.r.t. 
Kq can determine a reasonable order. So we have the follow- 
ing lemma. 



Lemma 2. A CNC is practically feasible if and only if the ET 
w.rt. Kq is acyclic. 



Proof: Let kd,eiz) = J2tLo kd,e,tz* and ye{z) = J^'^o Ue^tz^ 
be LEK for adjacent pair (c?, e) and the transmitted symbols 
over an outgoing channel e of a node v, respectively. Then 
form e G Out(y) by 



Veiz) ^^ye.tZ^ ^ X! ^dAz)yd{z). 
t=0 deln(v) 




Figure 5: The ET w.r.t. Kq in Fig. |3]is acyclic. 

Therefore, the transmitted symbol of e at time slot t is 

t 

ye,t = X! X! ^d.e,t-Tyd,T 
t 

= X! X! ^d,e,t-Tyd,T 

d£ln(v):d^e t=0 

+ ^ kd,e.t-Tyd,T 
deln{v):e^d r=0 

+ ^ kd,e,oyd,t- (2) 

de^/n(T;):e^d 

The data symbol ye.t can only be computed after the arrival 
of all yd,T, T < t and d £ In{v) at v. At time slot t, notice 
all yd,T, T < t have been received by v except those yd^t with 
can calculate y^ f according to (|2ji only if Kd,e.o — 
0. That is, a CNC is practically feasible w.r.t. {^i} only if 
kd,e,o = for all adjacent pairs {d, e) with e ^ d. 

Hence, according to the partial order {^,} defined on the 
acyclic ET, all yd.r, d < e have been calculated for all t < i 
because all incoming edges of v are previous to channel e. 
Meanwhile, all yd,T,s ^ d are also obtained from previous 
time slots t < t. That is, we can compute the output over 
channel e in terms of the first two items, whose elements are 
all known at time slot t. □ 

As for an acycUc network, the ET w.r.t. any LEK matrix K 
is acyclic. Moreover, K is a strictly upper triangular matrix, 
which is also a nilpotent matrix which is an nx n square matrix 
K such that iiT™ = for some positive integer matrix power 
m p2) . Meanwhile, it is known that the entry in if™ repre- 
sents the coding gain between any pair of channels through 
length m path. Hence there must exist an integer m < n such 
that if ™ = because the longest path is bounded in acyclic 
networks. However, over a cyclic network, the LEK matrix 
K{z) is no longer a strictly upper triangular matrix, we will 
discuss the properties of LEK in this case in the following sec- 
tion. 

III. The conditions to uniquely determine the 

GEK MATRIX F{z) 



So far, most work in convolutional network coding was de- 
veloped in the case of a unit time delay network or under the 
assumption that there is at least one delay along every cycle. 
In either case, the ET w.r.t. Kq of a given CNC is acycUc. We 
will show that Kq is nilpotent, which is also a necessary and 
sufficient condition to expand I + K{z) + K'^{z) + • • • . Then 
the GEKs can be uniquely determined by K{z) of the CNC 
based on Proposition 1 . 

Theorem 1. Kq is nilpotent if there is no cycle in the ET w.r.t. 

Kq. 

Proof: We prove by induction. Denote by if^" the entry of 
if™ in the i-th row j-th column, which indicates the coding 
factor between i and j along the paths of length m. The case 
for TO = 1 is trivial, where Kq represents the one hop trans- 
mission matrix. Assume that the statement holds for to — 1. 
By matrix multiplication, we can obtain 

I 

It is easy to see that at least one of K^i~^ and Ki j equals to 
zero for any I, because otherwise there must be a path from i 
to I of length (to — 1) and with I adjacent to j, which implies 
there is a path between i and j of length to. Given acyclic ET 
w.rt. ifo, the maximum paths between arbitrary two nodes are 
bounded. Hence A'o is nilpotent. □ 

The idea of this theorem is referred in |3!|. Inversely, a 
nilpotent Kq does not imply the acyclic ET w.rt. Kq. 




Figure 6: a cyclic network with if™ — 

Example 2. Assume that we multicast two messages from 
source S to both the nodes X and Y with given LEKs depicted 
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and Kq = 0. The corresponding ET in Figj?] contains two 
cycles. 




Figure 7: the ET w.r.t. K{z) in Figj6]is cyclic 

In the case of Figj6j F{z) is practically realizable, which 
is as a result of two overlapping cycles. The GEK for chan- 
nel 2 is derived from channels 1, 3 and 5. The same message 
components in channel 3 and 5 are counteracted. Therefore, 
the recursion along the cycle is destroyed, and the logical con- 
tradiction can be avoided. This is equivalent to the acyclic 
network coding in Fig|8] that is, channel 3 and 5 are actually 
not the incoming channels of node Y in the calculation of the 
output over channel 2. However, it remains unclear whether 
all the codes with nilpotent Kq can avoid the logical contra- 
diction. This will be left open for further work. 

We have shown that feasible F{z) yields to expandable 
/ + K{z) + K'^{z) + • • • . Now we will give an expandable 
condition of / + K{z) + /^^(z) H . 



equating the matrix coefficients, we obtain: 

Aq = /„ + Kq + kI + --- 

Ai = Ki + {KoKi + KiKo) 

+ {KoKoKi + KoK.Kq + K^KoKo) + ■■■ 
A2 = K2 + {KoK2 + K2Ko + KiKi) 

+ {KoKoK2 + KoKiKi + K0K2K0 + KiKqKi 

+KiKiKo + K2K0K0) + ■■■ 

JlH 



where Kj^, ^ Kq, 1 < jV < « is the coefficient matrix of 
ji' degree in z of K{z). It is obvious that ^0 is finite if and 
only if Kq is nilpotent. Let m be the smallest integer such that 
if™ = 0. Then mi' e {0, 1, • • • , m — 1} because the product 
term containing if™'' will vanish as mi' > m. To prove / + 
K{z) + K'^{z) + • • • is expandable, we have to prove Ai is the 
sum of finite terms for all i as if™ = 0. Now for every Ai, let 
us consider the number of the terms. Because ji + • • • + j; — i 
and 1 < ii' < i, the number of the solutions is less than i*, 
so there are less than choices of these Kj^, arrays. Then, 
for a certain Kj^, array, Kq can be divided by Kj^, into I + 1 
parts. So there are at most 771'+^ choices of Kq arrays because 
of mi' < TO, which is less than m*+^. Therefore, there are less 
than x m^'^^ terms for every Ai. □ 



Theorem 2. Let Kq be the constant coefficient matrix of 
K{z). Then I + K{z) + K'^{z) + • • • is expandable if and 
only if Kq is nilpotent. 

Proof: Assume A{z) - YZo = + E^i K"^^)- 




Figure 8: an equivalent encoder w.r.t. Fig|6] 
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Figure 9: realizable F (z) with Kq non-nilpotent 



The expandable / + K^z) + K'^(z) + • • • implies invert- 
ible /„ — K{z). But sometimes, /„ — K{z) may be invertible 
whereas Kq is non-nilpotent. 



Example 3. As in Figj9] the LEK matrix K{z) is 
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and for t = 0, Kq is 
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It is easy to check that Kq ^ 0, i.e., Kq is non-nilpotent. 

On the other hand, it is easy to know that /„ — Kq is invert- 
ible, and we have 



Fo = H,iI„-Ko)-' 

10 1111 
1110 



Meanwhile, F{z) can be derived from 

F{z) = 



1 1/(1 -z) 1 1/(1 -z) 1/(1 -z) 
1 1/(1 -z) 1 zfil-z) z/{l-z) 



We observe that the above Fq and the constant term in the 
matrix power series representation of F{z) are the same. 

Theorem 3. Consider an co-dim W-CNC with LEK matrix 
K{z). The fallowings are equivalent. 

1) In — Kq is invertible over F. 

2) In — K{z) is invertible over¥[{z)]. 

3) K{z) uniquely determines the GEK matrix F{z). 
Proof: 

1) => 2) Since /„ — Kq is invertible over F, then 
det (/„ — Kq) ^ 0. By using the following property of the 
determinant repeatedly. 
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&2 
022 



0-211 



a-nl 0,n2 



det {In-Ka-Kiz ) 

det (/„ -Ka)+z''Att{Ki) 



-z2"det {K2) 



+det (cross-term matrix) 

^ 

which indicates that there is at least one column of the cross- 
term matrix whose entries are divisible by z, therefore the de- 
terminant of the cross-term matrix is divisible by z. Because 
the constant term of det(/ — K{z)) is non-zero, / — K{z) is 
invertible over F[(z)]. 

2) =4> 3) Due to equation ([T|i, it is obvious that the GEK 
matrix F{z) can be uniquely determined by Hg{In—K[z))~^ . 

3) ^ 1) Let F{z) = Y.T Ftz' and K(z) = ^t^*- 
By substituting these into ([T]l, we have 



(3) 



Let the corresponding term coefficients be equal at both sides 
of the equation, then for t — 0, 



Fo{I-Ko)^H, 



(4) 



It is sufficient to show that F{z) cannot be uniquely deter- 
mined by K{z) if rank{I — Kq) < n. According to Fq is 
uniquely determined by K{z) only if its z-th row is the unique 
solution of Xi{I — Kq) — hi, i.e., rank{I — Kq) ~ n, where 
Xi and hi, i — 1, 2, • • • , w, are the i-th row of Fo and 
respectively. Therefore, / — Kq is invertible over F. 
Meanwhile, from ([3]), we obtain. 



Fi - FoKi - FiKo 
^ Fi 





FoKi{I-Ko) 



t-i 



Ftil ^ Ko) - } FrKt-r = 



T=0 



t-1 

Ft = {Y,FrKt-r){I-K„y 

T=0 



□ 

So far, we have discussed the conditions that the LEKs of a 
CNC can uniquely determine the GEKs. Fig{T]summarizes our 
results mentioned above. The conditions are characterized in 
terms of Kq. They simplify the convolutional network encoder 
design and provide a mathematical basis for CNC. 

IV. Decoding of CNC 
Consider a practically feasible w-dim F-CNC with GEK 



matrix F{z) = [/^(z)] 



at a sink node r. Similar to 



the field-based linear network codes, the F-CNC is decod- 
able at r if and only if rank{F{z)) = uj. That is, there are 



Lu incoming channels of r whose GEKs are linearly indepen- 
dent. Moreover, as adopted in |j5][TT][T8j, as long as we can 
find an x w decoding matrix D{z) over F[(z)] such 

that F{z)D{z) = z^I, the code is decodable with delay L. 
Based on D{z), the source symbol vectors can be sequentially 
decoded via the finite-state linear time-invariant shift regis- 
ters ||7j. However, the computation of D{z) or even the initial 
full rank check of F{z) over F[(z)] might encounter several 
issues in practice: 

• We need the full knowledge of F{z). However, due 
to possible existence of cycles in the network, even if 
the GEKs are polynomials, every entry in F{z) may 
involve infinite terms. When the CNC is deployed in 
a randomized manner, the information of F{z) will be 
carried along with the transmission symbols to the sink 
node time-slot by time-slot. Thus it will be impossi- 
ble to know all the terms in every entry of F(z) before 
checking its decodability, not saying the computation of 
D[z). 

• Given that F{z) is fully known to the receiver, the com- 
putation over F[(z;)] takes high computational complex- 
ities due to the possible occurrences of quotient of poly- 
nomials. 

• Given that F{z) is of full rank, we can not determine 
whether it is decodable with delay L untill the calcula- 
tion of D{z). 

The existence of a decoding matrix D{z) over F[(2:)] s.t. 
D{z)F{z) — z^I^ is a sufficient condition for the code's de- 
codability with delay L. In order to make the decoding process 
easier to be handled in practice, especially in the randomized 
settings, the main goal of this section is to formulate a series of 
necessary and sufficient conditions on the code's decodability 
with delay L, which will only deal with the first L -\- 1 terms 
Fo,...,FLmF{z)^Y.t>QFtzK 

Let x{z) = X^oo ^t-^* '■^^ power series of symbol vec- 
tors generated at the source and y{z) = X]t>o Vt^^ power 
series of symbol vectors received at sink node r. Moreover, 
denote by Xq{z) — X]f=o ^t-^* sequence of symbol vec- 
tors generated from time unit to L, and j/q (z) = J2t=a Ht^^ 
the sequence of received symbol vectors over the same period. 
Specific to the source node s and the sink r, a CNC on the 
network can be regarded as a linear sequential encoder with 
encoding kernel F{z). Thus, the definition of invertibility of 
linear sequential circuits in fl5] can be applied to physically 
define the decodability of a CNC. 

Definition 7. A practically feasible ¥-CNC is decodable at 
sink r with delay L if and only if for every non-negative inte- 
ger k, the input segment is uniquely determined by the 
response segment y^' 



■'^+\z). 



Theorem 4. A practically feasible ¥-CNC is decodable with 
delay L if and only if Xq{z) is uniquely determined by the re- 
sponse segment yQ{z). 

Proof Necessity follows directly from Definition |7] For suf- 
ficiency, suppose that Xq{z) = xq can be uniquely deter- 
mined by the sequence of output symbol vectors xoFq, xqFi + 

xiFo,xoF2-\-xiFi+X2Fo, ■ ■ ■ ,xqFl^ ^xlFq from time 

unit to L. At time unit L + 1, the output symbol vector is 
+ - • •+xl+2-Fo- After subtracting the linear gain of Xq 
from each of output symbol vectors from time unit 1 to L + 1, 
we can determine the source symbol vector xi in the same 
manner as xq from the sequence of modified output symbol 
vectors XiFq, xiFi + xjFo, • • • , XiF^ H h xl+iFq. Fol- 
lowing this way, at each time unit L+fc after subtracting the re- 
spective linear gains of xi, ■ ■ ■ , Xk~i from the output symbol 
vectors yL+i, ■ ■ ■ , VL+k, the source symbol vector Xk can be 
uniquely determined by the modified sequence of output sym- 
bol vectors XkFo, XkFi +Xk+iFo, • • • , Xfe^i H \-Xk+LFo. 

□ 

Since j/q (2) = '^^^q XiFL-iz\ whether the code is de- 
codable with delay L is only related to the first L + 1 terms in 
F{z). Thus, 

Corollary 1. A practically feasible CNC with GEK matrix 
F{z) at a sink r is decodable with delay L if and only if a 
practically feasible CNC with GEK matrix F^{z) at r is de- 
codable with delay L. 

The delay constraint L for the decodability of F^{z) in 
Corollary [1] is crucial. For instance, if the GEK matrix F{z) 
' 1 + z 1 + \ 

over F2, it is not decodable. However, 

1 1 + z J 
Fq{z) = ^ is decodable with delay 1. 

Justified by Theorem |4] as long as sink r is able to decode 
the first source symbol vector xq at a certain time unit L, it is 
able to decode the k*^ source symbol vector at time unit L + k, 
no matter how the received symbol vector yL+k is formed. For 
example, assume the GEK matrix at sink r is 



Fiz) 



z 

1 + z 



Since F{z) = Fo + F^z = (^^^ 1 ) + ( 1 ' 

easy to see that the symbol vector xq can be decoded at time 

1 



unit with the field-based decoding matrix Dq 



1 



Due to the linearity of CNC, it suffices to check whether 
Xq{z) is uniquely determined by 1/0(2;) in order to determine 
decodability with delay L. 



Assume that the input is x{z) = (1,0) + (1, l)z. Then the 
sequence of received symbol vectors is y{z) = x{z)F{z) = 
(1, 0) + (1, 0)z. At time unit 0, we can recover the first symbol 
vector Xq by 

yoDo = il,0) ( J ? ) =(1-0) 

At time unit 1, we know the linear gain of symbol vector xq in 
the received symbol vector yi is Fi = ^ ^ J ^ . The sub- 
traction of xqFi from j/i yields y[ = (1, 0) — (0, 1) = (1, 1). 



In the same way as getting xq, the second source symbol vec- 
tor xi can be recovered by 



(1,1)- 



This type of soft decision sequential decoding method has 
been adopted in Q . 

Now, the condition of decodability with delay L has been 
reduced to check whether xo{z) can be recovered from yg {z). 
Because 



viz) 



x{z)F{z) 



= ( 



E 



(5) 



we have y^(z) = Y.t=oiJ2i=o^i^t-i)z^- Alternatively, de- 
note by XL and y^, respectively, the row vectors {xq, ■ ■ ■ , xl) 
and (j/o, • • • tVl) of the data symbols generated at source s 
and received at r from time unit to L. Then, 



Vl = xlFl, 



where 



/ Fo 




F, 



Fi 
Fo 



V 



Fl 

Fl-1 

Fo J 



(6) 



(7) 



Based on Fl, which is a matrix over F, some necessary and 
sufficient conditions for recovering Xq{z) from (2) can be 
derived. 

Lemma 3. A practically feasible CNC with GEK matrix F{z) 
for a sink node r is decodable at r with delay L if and only 
if there is an \In{r)\ x w matrix D{z) over F[(z)] such that 

FlDl ~ ^ ^ ^ )' ^^^'^^ defined in the same man- 

ner as F Lfrom F{z). 

Proof: Since JJlDl = xlF^Dl ~ (0, • • ■ ,0,xq), the suffi- 
ciency follows. For necessity, if xq can be recovered from j/^, 

Dl 

then there is an {L + l)\In{r)\ x lj matrix D 



Do 



over F such that y^D — xo- Equivalently, 
( Fo Fl • • 



Fo 



V 



Fl 

Fl-1 

Fo ) 





( ) 






Dl-1 










) 


[ Do ) 





which implies that for every < j < i — 1, J2l>o FoF>i = 0. 
Therefore, for any x w matrix D{z) with first L + 1 

matrix terms Dq, ■ ■ ■ , in its matrix power series represen- 
tation, FlDl ( !! I . □ 





For technicity, assume F_i is a zero uj x |/n(r)| matrix. 



Theorem 5. A practically feasible uj-dim F-CNC with the 
GEK matrix F{z) for a sink r is decodable at r with delay 
L if and only if 



rank(FL) — rank{F l-i) = 



(8) 



Proof: We first show the necessity. According to Lemma [3] 
there is an x cj matrix D{z) over F[(z)] such that 

FlDl = 1 n I - Consequently, 







(9) 



and 



where D 



Dl 



{Fo ■■■ Fl)D = 
(0 Fl-i)D = {S. 

According to the matrix property 



(10) 



Do 



that rank[AB) < min{rank{A),rank{B)}, equation (j9]l 
implies rank{Fo ■ ■ ■ Fl) > cj- On the other hand, since 
(Fo • • • Fl) only has w rows, rank{Fo ■ ■ ■ Fl) < uj- We ob- 
tain rank{Fo ■ ■ ■ Fl) = cj. As a result, there does not exist a 
non-zero w-dimrow vector a subject to a- (Fq • • • Fl)D ^ 0. 
In other words, the intersection of the row space and the ker- 
nel of column space of D is only the zero vector. However, the 
row space of (0 F^-i) is a subspace of the kernel of the col- 



umn space of D because of ( 10 1. Therefore, there is no com- 
mon nonzero vector between the row space of (Fq Fi • • • Fl) 
and the row space of (0 Fl-i)- The necessity follows. 
For sufficiency, if rank{FL) — rank{F l-i) = then 

1 . the first lo row vectors in Fl have rank w; and 

2. there is no common nonzero vector between the vector 
spaces generated by the first a; and the last ujL rows in 
Fl- 

Denote by Null{0 Fl-i) the null space of last ljL rows in Fl 
and by (•) the row space of a matrix. Define a linear transfor- 
mation (f) : Null{0 Fl-i) with = (Fo • • • Fl)u. 
Prescribed by condition 2 above, 

{{Fo ••• Fl)){^{{QFl-i))^{0}. 
Moreover, the kernel of (f) is 

{w e F" : (0 Fl-i)u = and (Fq • • • Fl)u = 0}. 
Thus, 

dim(ker((/i)) 

= \In{r)\{L + l)-rank(FL-i)-rank{Fo ■ ■ ■ Fl) 
= \In{r)\{L + l) - rank(FL-i) - i^- 

On the other hand, 

dim(A?M//(0 Fl-i)) = \In{r)\{L + 1) - rank(FL-i)- 



Hence, the cardinality of the image of (j> is F"^ and then (p is 
surjective. As a result, there exist w \In{r)\{L + l)-dim col- 
umn vectors mi , • • • , over F such that 



(0 FL-i)iui 



□ 

By checking ([8]l iteratively, we can find the minimal decod- 
ing delay at sink r which has been characterized in p6) . If the 
CNC with GEK matrix F{z) for r is known to be decodable 
with delay L at r, an \In{r)\{L + 1) x u field-based decoding 
matrix as prescribed in Lemma [3] can be calculated from 
Fl. Based on this matrix, we can decode the source symbol 



vector Xk at time unit k + L,k > 0, as follows. 


Since 




f Fk ■■■ 


Fk+L \ 


(y/c, • • ■ ,2/fc+L) = (a;o, • • • ,Xk+L) 


Fo ••• 


Fl 




■■■ 






V ... 


Fo J 


{xk, ■ ■ ■ ,Xk+L)FL 








' Fk ... 


Fk+L \ 


= {Vk,--- -.yk+h) -xk-i j 








vFi ... 


Fl+L ) 



Thus, the source symbol vector Xk can be decoded via 

/ Fk ... Fk+L 



iVk 



,Vk+L) - Xk- 



\F, ... 
The field-based decoding algorithm ( 1 1 



Fl 



Dr 



decoding matrices subject to F^Dl 
instance, if the GEK matrix is F{z) — 



Do 
(11) 

may adopt different 






For 
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and Fl 



or 



, then either 









F2 
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^ J 
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1 
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1 





can be adopted to sequentially decode 



source symbol vectors xq, xi, ■ ■ ■ with delay 1. 

Similar to the one proposed in |10J , the present field-based 
algorithm is time variant. But the complete information on 
F{z) is not required here. On the other hand, if the matrix 
D{z) over F[(z)] is further computed such that F{z)D{z) = 
z^I^, the source symbol vector Xk can be time invariantly de- 



coded at time unit L - 
shift registers. 



k by, in particular, time invariant linear 



There is a useful consequence of Theorem |5] for prelimi- 
nary decodability check with lower computational complexity. 
This has been adopted in the adaptive random construction al- 
gorithm for a CNC in p7) . 

Corollary 2. A sink node r with GEK matrix Fr{z) = 
FfZ^ is decodable with delay L only if 



ik ( Fl 



Fl 



(12) 



This is a necessary but not sufficient condition to decode 
with delay L at sink r. When determining whether we can 
start decoding with delay L, we first check (12 1. From the 



time unit L that ([12]) is satisfied, we shall turn back to check 
dSll instead. 




Figure 10: A Convolutional Code 



Example 4. Fig 10 depicts the LEKs and GEKs of a practi- 
cally feasible CNC on a cyclic network. We shall illustrate the 
decoding procedure of this code at sink X. The GEK matrix 
for X is 



F{z) = 



1 z' 

z z' 

1 







1 



1 
1 



First, we check at each time unit i the rank of {Fo • • • 
until it is equal to w. An easy check in this example gets 
rank{Fo Fi) — uj. Then, starting from t = 1, we continue 
to check rank{Ft) — rank{Ft^i) until it is equal to lu. At 
t = 1, 



Fo 



/ 1 \ 

10 

10 

\ y 



Thus, rank{Fi) = 2 and rank{Fi) — rank{Fo) = 1. As a 
result, even if rank{FoFi) is of fuU rank, we are still not be 



able to start decoding. At time unit i = 2, we have 



Fq Fi F2 
Fo Fi 
Fo 



Hence, rank{F2) — 4 and 
implies that we are able to compute a matrix 



such that 
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000 
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rank[Fi) 


= 2. This 
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Following this, the source symbol vector xq can be decoded 
/ D2 

via (j/o ^12/'') Di I and Xk, k > 0, can be sequentially 

decoded via ([TT]|. 

On the other hand, if we plan to start the sequential de- 
coding xo,xi, - ■ ■ with delay 2 in a time-invariant manner, 
we need to find the matrix D{z) such that F{z)D{z) ~ z^I. 
Here, it can be deduced that 



D{z) = 



z^lil-z) -z^/{l-z) 
-z/{l-z) 1/(1 -z) 



V. Conclusion 



In this paper, we study convolutional network coding by 
means of matrix power series representation and conclude a 
few new results. Firstly, for a convolutional network code 
(CNC) over a single source network with possible cycles, we 
show that the nilpotent constant coefficient matrix of the 
LEK matrix K{z) is sufficient to determine whether the code 
is practically feasible. The encoding topology w.r.t. K^) is also 
introduced to illustrate the encoding order of the code. Addi- 
tionally, some equivalent conditions are presented to uniquely 
determine the GEKs based on K{z). 

For decoding of a CNC, we provide a physical definition 
of decodability at a sink node r with delay L, which only in- 
volves partial encoding information. Based on this new def- 
inition, several necessary and sufficient conditions are estab- 
lished for decodability according to the first L + 1 terms in the 
matrix power series representation of the GEK matrix F{z) at 
r. They yield a less computational method for decodability 
check, and does not require all terms of F{z), which may be 
infinite because of cycles in the network. As a result, CNC 
becomes possible to be deployed in a decentralized manner. 

There remain several open problems in the study of CNC. 
One of the major challenges is to design the coding scheme 



under the more practical scenario that delay conditions vary 
from time to time. The design of a CNC with minimal decod- 
ing delay is another interesting direction for future research. 
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